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Abstract 



O . 

■ Available Skyrme parametrizations with hyperons are examined from the point of view of 
their suitability for applications to neutron stars. It is shown that the hyperons can attenu- 
ate or even remove the problem of ferromagnetic instability common to (nearly) all Skyrme 
parametrizations of the nucleon-nucleon interaction. At high density the results are very sensi- 

^■f-^ ■ five to the choice of the A-A interaction. The selected parameter sets are then used to obtain 

I the resulting properties of both cold neutron stars and hot protoneutron stars. The general 

features known from other models are recovered. 

> '. PACS numbers: 97.60.Jd, 21.65. +f, 21.30.-x 

■ keywords: neutron stars, Skyrme forces, hyperons 

o 

O ' 1 Introduction 
^ ■ 

I Skyrme parametrizations provide a simple tool for calculating the equation of state. Skyrme or 
Skyrme-inspired models are often used in applications to neutron stars with adapted versions which 
pay special attention to the asymmetry properties. The Skyrme Lyon series [1] is the most modern 
O I example of this approach and was designed to reproduce saturation properties, the main properties 
of nuclei as well as the results of simulations of pure neutron matter. 

In neutron stars in chemical equilibrium, hyperons appear at a threshold density about twice 
the nuclear saturation density. A sizable hyperon fraction can be expected to be present in the 
^ ' core of the most massive neutron stars; the central density of the Skyrme Lyon series for example 

■ is of the order of seven times saturation density in the star with maximum mass and three to four 
times saturation density in a 1.4 Mq star. 

Original Skyrme parametrizations did not take hyperons into account. Balberg and Gal [2] have 
remedied the fact with their parametrization of the energy density. It could be argued however 
that the description of the nucleon sector provided by this parametrization is rather poor. 

Rather than designing a new force in the hyperonic sector with the purpose of describing neutron 
stars, an other possibility would be to restrict oneself to the numerous available data sets which 
were fitted to nuclear properties. This is the approach taken by Rikovska Stone et al. [3], who 
recently tested 87 parameter sets for the NN force in neutron star calculations. 

In the hyperonic sector there exist a few parametrizations of the Lambda-Nucleon and Lambda- 
Lambda interaction fitted to the properties of hypernuclei, including several well tested sets by 
Lanskoy el al. [4, 5] and an earlier set by Fernandez et al. [6]. The present work can be considered 
as an extension of the study of Rikovska Stone et al. including the A hyperon. 
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A general feature of Skyrme models in the np sector is that they show a ferromagnetic transition 
at rather low density [7, 8]. If present, such a transition, apart from modifying the equation of 
state and chemical equilibrium, could give rise to induced magnetic fields in rotating neutron 
stars. It also plays a crucial role in calculations of the neutrino mean free path in supernova and 
protoneutron stars. Calculations of this parameter in non relativistic (Skyrme or Gogny) models 
of nuclear matter [9, 10, 11] thus obtained that the cross section diverges and the mean free path 
drops to zero at the transition. 

It remains an open issue whether this transition is genuine or an artefact of the Skyrme model. 
While a free Fermi gas eventually becomes ferromagnetic, the nuclear correlations are known to 
play a crucial role. Most non relativistic calculations which take correlations into account, e.g. by 
solving the Brueckner-Hartree-Fock equations with modern bare NN potentials [12, 13], concluded 
that spin ordered matter was not favored energetically. Relativistic mean field models predict no 
ferromagnetic transition; on the other hand relativistic calculations in the Hartrec-Fock approxima- 
tion by Bcrnardos et al. [14] find a transition, albeit as a rather large density (4 rigat)- A recent work 
by Maruyama and Tatsumi [15], although not putting forward a quantitative prediction, reaches a 
similar conclusion as to the importance of the Fock contribution. 

In this work we will take the viewpoint that the ferromagnetic transition, especially when 
appearing at such low densities as 2-3 times saturation, is probably an artefact of the Skyrme 
model, and will select among the available Skyrme parametrizations those which as far as possible 
avoid or delay the transition to higher densities. 

This work includes the effects of temperature and neutrino trapping and can also be applied in 
hot protoneutron stars. 

A further motivation of the present study was to obtain a model which would allow to calculate 
the neutrino-baryon scattering rate in the hot protoneutron star formed shortly after the supernova 
collapse. It was therefore necessary to be able to calculate the Landau parameters in the spin 5 = 1 
channel. Besides determining the position of an eventual ferromagnetic transition, the Landau 
parameters in the spin S = 1 channel play a central role when calculating the axial response 
function in the random phase approximation, which is the dominant contribution to the neutrino- 
baryon scattering rate. This application is presented in a separate paper [16]. 

After a presentation of the available Skyrme parametrizations in section §2, the following sec- 
tions proceed to select a few among them according to their suitability for neutron star applications. 
Section §3 first examines the threshold density for hyperon formation under the conditions of (3 
equilibrium. §4 presents a discussion of the issue of fcrromagnctism. The selected parameter sets 
are used in section §5 to obtain the properties of the corresponding neutron stars. The effect of 
non vanishing temperature and of a trapped neutrino fraction are studied in §6. The main results 
are collected and discussed in the conclusion. 

The formalism is kept to a minimum in the main text, while all relevant formulaes are gathered 
in the Appendix. 

2 Skyrme parametrizations of the hyperonic sector 

The usual Skyrme model of nuclear matter introduces the nucleon-nucleon potential 

VNNiri - ra) = to (1 + xoP„)5{ri - ra) + ^ti (1 + xiP„) [k'"^ 5{ri - rs) + <5(ri - rg) k' 

+t2 (1 + X2Pa)k' 5{ri -r2)k + ^t^ (1 + x^P^) p% (^-^^^ <^(n - rs) (1) 
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It can readily be generalized to include nucleon-Lambda and Lambda-Lambda interaction potentials 
(see e.g. [4, 5]) 



VNA{rN - va) 



uo (1 + yoPa)S{rN - ta) + -^1 fc'^ S{rN - ta) + d{rN - va) k"^ 




VAK{ri - r2) 




The potentials normally also include spin orbit contributions. They are not explicited here since 
they will not be used in the remainder of this paper. 

The energy density is then obtained in the Hartree-Fock approximation from 



In homogeneous matter the wave functions are formed from antisymmetrized plane wave states. 
The explicit expression of the energy density is given in the Appendix. 

In our search for some suitable sets of Skyrme parametrizations we have tested 43 NN Skyrme 
forces in combination with 13 parametrizations of the NA and 4 options for the AA forces. The 
NN forces were chosen among: SIII, SkM*, SI', SIIF, SV, RATP, SGI, SGII, SLy230a, SLy2, SLy4, 
SLy5, SLy6, SLy7, SLy9, SLylO, SkO, SkO', SkS3, Skll, SkI2, SkI3, SkM, SkI5, Rs, Gs, SkT4, 
SkT5, T6, SkP, Ska, MSka, SK272, SK255, Skz-0, Skz-1, Skz-2, Skz-3, Skz-4, SkSc4, SkScl5, MSk7, 
SKX. After demanding that these sets fulfill various constraints that will be discussed in detail in 
the next two sections, only four NN forces were sorted out: the SLylO force from the Skyrme-Lyon 
series [1], the modern SkI3 and SkI5 forces by Reinhard and Flocard [17] and the older SV force 
by Bcincr et al. [18]. 

For the nuclcon-A interaction we have tried several parameter sets given by Lanskoy and Ya- 
mamoto [4, 19] and by Fernandez et al. [6]: The A'^A sets numbered from I to V in [4] are named 
LYI-LYV here, the sets numbered from 1 to 6 in [19] are named YBZ1-YBZ6, other choices are 
the SkSHl and SkSH2 sets of [6] or to switch off the NK interactions. Finally, Lanskoy gives in [5] 
three sets SLLl, SLL2, SLL3 for the AA interaction, or we can switch it off. 

The parametrization of Lanskoy and Yamamoto [4] was extracted from G-matrix calculations 
[20] performed with the Jxilich and Nijmegen potentials and were tested on hypernuclei. Other NK 
and AA interactions were fitted directly to hypernuclei data. The values given by Lanskoy and 
Yamamoto [4] assume that the nucleon sector is parametrized by the SkM* or SIII interactions 
while the parameter set of the Salamanca group [6] was used together with the SkS3 interaction. 
Even though this is not fully consistent, we also used other more modern parametrizations in the 
nucleon sector like the Skyrme Lyon sets [1] in order to investigate the role of the ferromagnetic 
transition (see section §4) 

Selected values of the parameter sets are listed in tables I, II and III. For other parametrizations 
we refer the reader to the exhaustive table published by Rikovska Stone et al. in the case of the 
NN forces and the papers of Lanskoy et al. and Fernandez et al. for the NA and A A forces. 

As this model only includes the A hyperon, it appears to be less complete than the model of 
Balberg and Gal [2] which includes all the hyperons. As mentioned in the introduction, we fixed 
our choice on the parametrization by Lanskoy et al. because these authors provide the two-particle 
potential rather than only the unpolarized energy density. This leaves open the possibility of going 




= £nN + £nA + ^AA 



(4) 
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beyond the Hartree-Fock approximation and to investigate the response functions of the matter at 
the RPA level. 

It would be straightforward to extend the model to take into account other hyperons such as 
the Data about the S~ in nuclear matter however is rather scarce and the author is not 

aware of an equivalent Skyrme parametrization including the S at the same level of precision. 
One possibility would be to adjust Skyrme parameters to reproduce Brueckner calculations e.g. 
using density matrix expansion techniques. On the other hand the neglect of the S hyperon is 
sometimes justified from the very lack of observation of E hypernuclei, when this fact is interpreted 
as indicating that the E-nuclcon force is in fact repulsive. In that case /5-equilibrium equations 
would predict that the threshold for E hyperon formation in neutron stars is shifted to very large 
densities (see e.g. [2, 21, 22, 23]). 

We will use the model of Balberg and Gal [2] and another by Banik and Bandhyopadhyay [24] 
in order to estimate the error committed by neglecting other hyperons and in particular the E~. 



Model 


Q 


/() 


/i 


l-i 


I 'A 


,ro 


,ri 


•T2 


J":-! 


SLylO 


1/6 


2506.77 


430.98 


-304.95 


13826.41 


1.0398 


-0.6745 


-1.0 


1.6833 


SkI3 


1/4 


-1762.88 


561.608 


-227.09 


8106.2 


0.3083 


-1.1722 


-1.0907 


1.2926 


SkI5 


1/4 


-1772.91 


550.84 


-126.685 


8206.25 


-0.1171 


-1.3088 


-1.0487 


0.3410 


SV 


1 


-1248.3 


970.6 


107.2 


0. 


-0.17 


0. 


0. 


1. 



Table I: Skyrme parameters for the NN interaction 

(to is given in MeV.fm^, ti and t2 in MeV.fm^, in MeV.fm^+^"^, the other parameters are 
adimensional) 



Model 




no 


Ul 


U2 


W3 


2/0 


2/3 


Vk 


LY-I 


1/3 


-476. 


42. 


23. 


1514.1 


-0.0452 


-0.280 


-27.62 


LY-IV 


1/3 


-542.5 


56.0 


8.0 


1387.9 


-0.1534 


0.1074 


-28.17 


YBZ-1 


1 


-349.0 


67.61 


37.39 


2000. 


-0.108 


0. 


-26.52 


YBZ-5 


1 


-315.3 


23.14 


-23.14 


2000 


-0.109 


0. 


-28.50 


YBZ-6 


1 


-372.2 


100.4 


79.60 


2000. 


-0.107 


0. 


-24.98 


SkSHl 




-176.5 


-35.8 


44.1 










-27.68 



Table II: Skyrme parameters for the NA interaction [4, 19] 

{uq is given in MeV.fm^, ui and U2 in MeV.fm^, -^3 in MeV.fm^"*"^^, the other parameters are 
adimensional. V\ is the potential felt by a A hyperon in nuclear matter at saturation and is given 
in MeV.) 



Model 


Ao 


Ai 


SLLl 
SLL2 
SLL3 


-312.6 
-437.7 
-831.8 


57.5 
240.7 
922.9 



Table III: Skyrme parameters for the AA interaction [5] 
(Ao is given in MeV.fm^ and Ai in MeV.fm^) 

3 Threshold for hyperon formation 

In neutron star matter we impose that the conditions for /^-equilibrium are fulfilled. The neutron, 
proton and Lambda hyperons are subject to the processes 

p + e" ^ n + Ve ; p + e" ^ A. + Ue 
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as well as their inverse 

n ^ p + e~ +Ve ; A ^ p + e'' +Ve (5) 
We write the equality of the chemical potentials 

A = (/Un + mn) - {jlp + nip) = fie- IJ'P ; + irin = IJ,A + mA (6) 

We must also impose electric charge conservation 

rie = rip (7) 

At a given baryonic density the electron, proton and hyperon fractions arc determined by the 
solution of Eqs. (6,7). The chemical potentials are determined by deriving the Skyrme energy 
density functional with respect to the density of corresponding particle. Their explicit expression 
is given in the Appendix. The electrons are relativistic and their chemical potential is given by 
He = \l + rrig. In protoneutron star matter with trapped neutrinos we have fii, = (67r^ni/)^^/^^, 
while in colder neutrino- free neutron star matter /ij^ = 0. 

The equations should actually take into account the muons in the equation for charge conser- 
vation together with the condition fj,e = IJ-n- The muons appear namely around saturation density. 
The muons were neglected in this simplified model. Their effect is in fact not very important, 
especially in our case where the charged S~ is also absent from the model. 

The threshold for A hyperon formation is determined by the condition 

HA{thr) = HAikpA = 0) = /u„ -h m„ - ttt-a 

= (l + f) Piv + (l + f ) + I [ni{2 + yi) + ^2(2 + ^2)] kp%' (8) 

with k = (3/5) (Stt^)*-^/^-*. It depends on the strength of the AN force through the parameters 
uo,ui,U2,U3.. and indirectly on the NN force through the chemical potential of the neutron (taken 
in npe matter in /3 equilibrium) and is independent of the A-A force. Depending on the choice of the 
parameters it is found that the threshold generally occurs between 1.7 risat and 4 nsat in agreement 
with other Brueckner-Hartree-Fock or relativistic mean field calculations. When numerical values 
are inserted one can see that the value of ^Aithr) is the result of a delicate cancellation between 
the no and ^3 terms, the contribution of the ui,U2 term being of the order of the sum of the uo and 
W3 terms. For example at pN = 3 ngat and rigat = 0.16 fm~^ we have for the LY-I parametrization 
of the NA force 

uo (1 + y^'^pN = -223.32 MeV , ^ua (l + ^) = 183.51 MeV 

I [ui{2 + yi)+ U2{2 + 2/2)] kp^"^ = 27 MeV 
o 

Let us also note that the value of pA^kpA = 0) as saturation density with equal number of 
neutron and protons is nothing but the single particle potential felt by a A impurity in nuclear 
matter, i.e. it should equal to the binding potential Va — —28 MeV obtained from data on 
hypernuclci. The actual values of = p.A{kFA = 0,nB = 0.16 fm~^) are reported in the last 
column of table II. 

We may conclude this section by stating that, once the parametrization of the NA force has 
been chosen so that the potential felt by a A in nuclear matter reproduces the experimental value, 
the condition that the threshold density for A hyperons in nuclear matter in /3-equilibrium should 
lay around 2-3 times saturation density does not severely constrain the admissible Skyrme NN 
and NA forces. We observe that stiffer equations of state (eos) and the eos allowing for a larger 
proton fraction have lower hyperon thresholds (see table V for a sample of the results) . 
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4 Transition to ferromagnetic state 



4.1 Criterion for a ferromagnetic instability 

Previous studies on the neutrino mean free path in neutron matter [9] or npe" matter in /3 equi- 
hbrium [11] found that a pole appears in the calculation of the axial structure function above a 
certain critical density. This feature is typical of Skyrme models and is related to a transition 
to a ferromagnetic state. In this section we will study how this critical density is affected by the 
presence of hyperons. 

Let us define the magnetic susceptibilities Xij where {n,p,A}: 

^rOM^r ^^ = ^^^^^T-m) (9) 

where £ = £{Pnlj Pnl^ Ppt^ Ppl^ PAi, PAi) is the polarized energy density functional, Aii are the mag- 
netizations and Pi are the magnetic moments. The inverse susceptibilities are therefore proportional 
to the second derivatives of the energy density functional with respect to the polarizations: 

l^il^j 2p ld'^{£/p) Pi]- Pa 

Xij PiPj 2 dsidsj pi^ + pii 

On the other hand it can be shown that the Ajj are related to the Landau parameters gQ defined 
in the Appendix through 



Pi -'^0 

A criterion for the appearance of the ferromagnetic phase is that the determinant of the inverse 
susceptibility matrix vanishes 

/l/Xnn l/Xnp l/XnA\ 

Det 1/xpn l/Xpp l/XpA = (12) 

\l/XAn 1/XAp l/XAA/ 

and in terms of the Landau parameters: 



Det 



/ (1 + Gl^^) 

(l + G^ 1=0 (13) 




It can be shown that this quantity appears in the denominator of the static axial response (see 
[16]). 



4.2 Dependence of the criterion on pEirameter sets 

A few NN forces (MSk7,SkX,SkS3) had to be discarded as not convenient for obtaining the equation 
of state away from saturation, some more for not permitting the formation of hyperons below 5 risat, 
others again for undergoing a transition to pure neutron matter before the threshold for hyperon 
formation was reached (SIII, SI', SIII', Skz0-Skz4). 
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The remaining ones all display a transition to the ferromagnetic state at some critical density 
n^""^*^ when no hyperons are present, with the exception of the SV force. Usually the ferromagnetic 
transition occurs earlier in pure neutron matter (PNM) than in symmetric nuclear matter (SNM), 
and at an intermediate density for npe matter in /? equilibrium. One exception is the case of the 
Lyon Skyrme forces where the reverse situation appears to occur. A closer look nevertheless reveals 
that only traces of protons are enough to drastically lower the critical density in very neutron rich 
matter so that the usual pattern is in fact recovered. 

Some forces (SkO, SkO', Skll, SkI4, SkT4, SkT5) again need to be discarded since the ferro- 
magnetic transition occurs below nsat in symmetric nuclear matter, so that the nuclei would in 
fact be unstable with respect to spin fluctuations. This leaves us with forces for which the n^^^^ 
occurs in the range ngat ~ 3 ngat) and the threshold for A hyperon formation in the range 1.7 ngat ~ 
4 nsat- A key point is now whether n^~^^^ is greater or lower than n^^. While the criterion for the 
onset of the ferromagnetic instability Eq. (13) decreases with increasing density in npe matter in 
(3 equilibrium, we noted that it always tends to increase again when enough hyperons are present. 
This then rejects the critical density n^^^^^^ for the onset of ferromagnetism in npAe matter in /3 
equilibrium to higher densities. In a few cases the instability even disappears altogether. For the 
majority of NN parameter sets, the ferromagnetic transition occurs before the threshold for hy- 
peron formation. In some cases (SLy4, SLy7, SGI) we have n^j. ~< n^rrT^' ^^"^ ferromagnetic 
transition is only delayed by a tiny amount. 

We finally selected four NN forces (SLylO, SkI3, SkI5, SV) as relevant for our purpose, namely 
to study neutron star matter with an hyperonic component before the ferromagnetic transition sets 
in. In Figure 1 (a c) wc have represented the criterion (13) for these parametrizations of the NN 
Skyrme force and various choices of the NA and A A forces. The modification of the slope at the 
hyperon threshold around 2 Ugat is clearly visible on Fig. 1-a. 

The SV NN interaction is atypical as it has no ts term. This parametrization is especially 
interesting since, among all tested NN interactions, it was the only one which does not give a 
ferromagnetic instability in npe matter in /3 equilibrium. This parametrization is rather old (1975) 
but yields an acceptable description of the properties of nuclear matter and nuclei and it complies 
with all the conditions necessary for describing a viable neutron star. The SLylO, SkI3 and SkI5 
are modern forces. The SLylO set was designed to reproduce features of pure neutron matter as 
obtained from the variational calculations of Wiringa et al. in view of its application to neutron 
stars. The SkI3 and SkI5 sets were designed in order to improve isotope shifts and incorporate the 
dipole sum rule enhancement factor k = 0.25. Let us notice that the SkI3, SkI5 and SV allow for 
larger proton fractions than the SLylO, in particular we can see from table IV that the criterion 
for opening of the direct URCA process with nucleons is reached before the theshold for hyperon 
production. 

The behaviour of the criterion for ferromagnetic instability is qualitatively very similar for all 
choices of the NA interaction. Forces yielding smaller pA{thr) and softer eos (see §5) also tend to 
be more efficient in lifting the ferromagnetic criterion det[l + Gij] above the critical zero axis. This 
is especially the case for the YBZA parameter set which permitted to avoid the pole also in the 
SLy4 and SLy7 parametrizations of the NN force. However the YBZ4 parameter set was rejected 
by Lanskoy on the ground that it gives overbinding of the A in hypernuclei so that we will not 
consider it further. Despite this overall similarity the actual presence and position of the pole is 
sensitive to the choice of the NA force. This happens because the det[l + Gjj] has already decreased 
considerably prior to p\{thr ) and the turnover due to the contribution of the A therefore must take 
place in the vicinity of the zero axis. 

We can also see on this figure that the A — A interaction plays an important role in determining 
^fen-T*^^- '^'^^ SLL3 (which also gives a stiffer eos, see next section) is less efficient in removing 
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npeA, SLL3 • 



Fig. 1 - Criterion for the transition to a 
ferromagnetic state. 



the pole while the SLLl set which gives the softest eos is also the most efficient in preventing the 
onset of ferromagnetism. The A- A force of Lanskoy [5] is very schematic (it is parametrized by 
Ao and Ai only, see the Appendix); moreover it was adjusted to the older value Ai?AA = —4.8 
MeV instead of the value AB\\ ~ — 1 MeV recently extracted from the "Nagara event". The case 
where a vanishing A-A interaction is assumed may in fact be closer to the true situation. A better 
knowledge of the A-A interaction at high density is needed. 

The values of density and proton content at the threshold for hyperon formation as well as the 
critical values for the ferromagnetic transition for various models are gathered in tables IV a-d. All 
densities are quoted in units of the saturation density of the corresponding model. The mention 
"grazing" means that det[l — Gij], while not actually crossing the zero axis, comes so near it that 
for all practical purposes the static axial response function will behave as if a pole were present. 

a) NN force=SLylO 

[threshold for ferromagnetism in PNM: 3.837 Wsat; in SNM: 3.821 ngatj in npe matter in 
/^-equilibrium: 3.055 risat] 



NA 


«thr 




no AA 


SLLl 


SLL2 


SLL3 


LY-I 


2.719 


0.041 


grazing [4.3] 


no pole 


6.183 


3.898 


LY-II 


2.868 


0.040 


3.160 


3.212 


3.236 


3.328 


LY-IV 


2.764 


0.041 


3.657 


no pole 


6.480 


3.818 


YBZ5 


3.319 


0.037 










YBZ6 


5.101 


0.029 










SKSHl 


2.294 


0.044 


4.577 


5.297 


4.729 


4.153 
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b) NN force=SkI3 

[threshold for ferromagnetism in PNM: 2.298 rXgat) in SNM: 5.73 ngat; in npe matter in 
/3-equiHbrium: 3.078 risat] 



NA 


' ''tlir 


Y 


no AA 


STJ,1 

kJ J-J J-J -L 


SLL2 


STJj3 


LY-I 


1.864 


0.126 


no pole 


no pole 


9.234 


4.663 


LY-IV 


1.873 


0.127 


no pole 


no pole 


8.900 


4.548 


YBZl 


2.015 


0.139 


no pole 


no pole 


8.753 


4.893 


YBZ5 


1.942 


0.133 


no pole 


no pole 


7.550 


4.080 


YBZ6 


2.076 


0.144 


no pole 


no pole 


9.194 


5.633 


SKSHl 


1.739 


0.116 


8.225 


9.495 


6.765 


4.425 


SKSH2 


1.681 


0.111 


18.875 


15.250 


8.253 


4.212 



c) NN force=SkI5 

[threshold for ferromagnetism in PNM: 1.772 ngatj in SNM: 2.659 ngatj in npe matter in 
/3-equilibrium: 2.140 risat] 



NA 


^thr 




no AA 


SLLl 


SLL2 


SLL3 


LY-I 


1.727 


0.150 


grazing [4.34] 


no pole 


4.707 


3.166 


LY-IV 


1.734 


0.150 


grazing [4.14] 


no pole 


4.627 


3.095 


YBZl 


1.838 


0.162 


no pole 


no pole 


8.656 


3.622 


YBZ5 


1.785 


0.156 


no pole 


no pole 


5.404 


3.239 


YBZ6 


1.882 


0.167 


no pole 


no pole 


9.118 


3.950 


SKSHl 


1.629 


0.138 


3.337 


3.914 


3.541 


3.103 



d) NN force=SV 

[threshold for ferromagnetism in PNM: 4.850 rigat) in SNM: no pole, in npe matter in 
/3-equilibrium: no pole] 



NA 




Yp 


no AA 


SLLl 


SLL2 


SLL3 


LY-I 


1.793 


0.125 


no pole 


no pole 


9.330 


4.576 


LY-IV 


1.800 


0.126 


no pole 


no pole 


9.143 


4.533 


YBZl 


1.906 


0.133 


no pole 


no pole 


8.649 


4.635 


YBZ3 


1.672 


0.116 


no pole 


no pole 


11.700 


5.050 


YBZ5 


1.851 


0.129 


no pole 


no pole 


7.521 


3.956 


YBZ6 


1.951 


0.136 


no pole 


no pole 


9.093 


5.200 


SKSHl 


1.691 


0.118 


no pole 


no pole 


8.920 


4.250 


SKSH2 


1.635 


0.114 


no pole 


no pole 


8.253 


4.068 



Tables IV a— d: thresholds for ferromagnetism in npAe matter in /3 equilibrium 

5 Equation of state and neutron star structure 

5.1 Equation of state, effective masses and the hyperon fraction 

For neutron star matter in beta equilibrium without hyperons, the equations of state are in or- 
der of decreasing stiffness parametrized by the NN interaction SV > SkI3, SkI5 > SLylO. The 
parametrizations SkI3 and SkI5 yield nearly indistinguishable results, so that SkI5 will not con- 
sidered further in the remainder of this work. When the hyperons are taken into account in the 
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calculation of the (5 equilibrium, the equation of state softens as expected. For a given NN interac- 
tion, we have in order of decreasing stiffness YBZ6 > YBZl > LY-I, LY-IV > YBZ3 > SKSH2 > 
SKSHl. Finally, when varying the AA interaction for given NN and NA forces, we obtain SLL3 > 
SLL2 > SLLl. If the AA interaction is set to zero, the equation of state is somewhat stiffer than 
SLL2 at ns < 5 rigat and much softer above 5 ngat- We note the sizable effect of the A- A interaction 
on the equation of state, as well as in the previous section on the value of the Landau parameters. 




The equation of state is shown on Fig. 2 for two of our preferred NN forces, and with the same 
choice LY-I for the NA interaction and SLL2 for the AA interaction. Even though SkI3 was stiffer 
that SLylO without hyperons, the combination SkI3-|-LYI-|-SLL2 gives rise to a higher hyperon 
content at a given density (see also Fig. 4) than SLylO-|-LYI-|-SLL2, so that it is also subject to 
more softening. As a consequence, the eos with hyperons are very similar. 

Figure 3 illustrates the behavior of the effective masses in the case of the parameter set SkI3 
-I-YBZ6-I-SLL2. The effective mass of the neutrons are lower than that of protons in neutron rich 
matter, a feature generally encountered in Skyrme models. The presence of hyperons cause both 
the neutron and proton mass to decrease less rapidly at high density. 




Figure 4 shows the particle fractions 1^ = rii/nB as a function of total baryonic density ub 
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for the SkI3+YBZ6+SLL2 and SLylO+LYI+SLL2 parameter sets. It can be seen on the left 
panel that the hyperons appear at higher density and are less numerous for a given density with 
SLylO+LYI+SLL2 than with SkI3+YBZ6+SLL2. The right panel shows the chemical composition 
of matter with a non vanishing number of trapped neutrinos at zero temperature. The case with 
trapped neutrinos and finite temperature relevant for protoneutron stars will be discussed further 
in §6 

5.2 Solution of the Tolman-Oppenheimer-Volkoff equation 

Our selected parameter sets still have to pass the test of causality in the density range of interest, 
and whether they can support neutron stars with maximum masses larger than the observed value 
1.4 Mq. 

The properties of neutron stars formed of npe matter in /? equilibium were calculated by 
Rikovska Stone et al [3] for a large number of Skyrme parametrizations. The four NN interac- 
tions that we singled out in section §4 all belong to the subset of interactions selected by these 
authors as giving viable neutron stars, with maximum masses of the order 2 Mq to 2.4 Mq. The 
central density reached in stars with the maximum mass is slightly larger than the density at which 
the velocity of sound reaches the velocity of light for the SLylO and SV parametrizations whereas 
the equation of state obtained with SkI3 and SkI5 always remains causal. In any case, stars with 
the fiducial mass 1.4 Mq always fulfill nc(1.4 M^^ < nsic^ = 1). 

When hyperons are added, the equation of state being softer, the limit = 1 is reached for 
larger densities. It is pushed to ~ 8 - 11 times saturation density for the SLL3 parametrization of 
the A-A interaction, to ~ 13 - 16 rigat for the SLL2 parametrization and is larger than 20 ngat for 
the SLLl set. The case where a vanishing A-A interaction is assumed always remains causal. In any 
case, such densities are beyond the range of validity of the model. As a consequence, although the 
softening of the eos with hyperons gives rise to higher compression rates and larger central densities 
in neutron stars than in stars made of npe matter only, the criterion < 1 is always fulfilled up 
to the central density of the most massive stars in our neutron star models with hyperons, except 
a few instances involving the SLL3 parameter set. 

We have solved the Tolman-Oppenheimer-Volkoff equation to obtain the (non-rotating) neutron 
star masss-radius relation. The equation of state is matched at lower densities with that of Negele 
and Vautherin [25] for p G [0.001 - 0.08] fm"^ and with the Baym-Pethick-Sutherland (BPS) [26] 
equation of state for p < 0.001 fm~^. The results^ are summarized in Table V a-d and a sample 
of mass-radius curves is shown in Fig. 5 and 6. For each combination of the NN, NA and AA 
interaction we give the density at which the speed of sound becomes superluminal, the central 
density and radius of a 1.4 Mq in case it can be formed, and the central density and radius of the 
star with maximum mass Mmax- 

As a rule of thumb it is known that the stiffer the equation of state is, the larger will be the 
maximum mass of the neutron star. Our results for Mmax follow accordingly to the classification 
in stiffness given at the beginning of this section. Thus the SKSHl choice generally yields Mmax < 
1.4 Mq whereas the YBZ6 always pass this test successfully. We see again that the choice of the 
A-A interaction is crucial to determine whether or not the star is able to reach to 1.4 Mq line. 
When the A-A interaction is switched off or for the SLLl choice, the equation of state is generally 
too soft to support a 1.4 neutron star. On the other hand the SLL3 always succeed in producing 
-^max > 1.4 Mq. The SLL2 choice barely makes it to 1.4 Mq for intermediate choices of the NA 
interaction. This results in a plateau feature, where tiny additions of mass effect a large reduction 
of radius without actually reaching the instability, the M = 1.4 Mq point being finally reached at 

^The parameters we obtained for the stars composed of npe matter may differ very slightly from those quoted by 
Rikovska Stone et al., presumably due to a different matching at low density or the neglect of the muons. 
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Fig. 5 - Neutron star mass-radius rela- 
tion with the Skyrme SkI3-|-YBZ6 inter- 
action, with and without hyperons. 
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unprobably high densities p > 10 ngat much beyond the vahdity of the model. A stiff equation of 
state on the other hand is not convenient since it favours ferromagnetism to appear earlier. The 
latter condition usually rules out the SLL3 choice. 

For the SLylO choice of the NN interaction ferromagnetism always sets in before the central 
density of a 1.4 Mq star is reached. In particular, it must be reminded at this point that ferro- 
magnetism is reached before the threshold of production of hyperons for the choices YBZl, YBZ5, 
YBZ6. For the SkI3 choice, the criterion nc(1.4 Mq) < riferro is satisfied for the sets (SkIS -|- LY-IV 
+ SLL3), (Skis + YBZ5 + SLL3) and (SkIS + YBZl or YBZ6 + any choice of AA ). The star with 
the maximum mass also fulfills all requirements if it is described by (SkIS + YBZ6 + no AA ) or 
(Skis + YBZ6 + SLL2). The results obtained with the choice SkIS are very similar. The choice SV 
represents an intermediate situation. The 1.4 Mq star fulfills both criteria nc(1.4 Mq) < nB{cl = 1) 
and nc(1.4 Mq) < nferro in the cases (SV + LY-I + SLLS), (SV + LY-IV + SLL2), (SV + YBZl 
+ SLL2/3), (SV + YBZ6 + any choice of AA). Moreover the star with the maximum mass also 
fulfills both criteria for the choices (SV + YBZ6 + no AA) and (SV + YBZ6 + SLL2). 
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a) NN force=SLylO 

[in npe matter in /3-cquilibrium: ^^(c^ = 1) = 7.308 n^at, ^^{lA Mq) = 3.60 ngat; 
R{1A M0)=11.O5 km, n,„ax = 7.69 nsat, M„iax = 1-99 Mq, i?(Mmax) = 9.52 km] 
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b) NN force=SkI3 

[in Ttpc ina^ttGr in /3-Gquilibriuni: '^sis^ ~ ~ 6.343 Tt^oXt 

(1.4 Mq) = 2.27 nsat, 

R{IA Mq)=U.21 km, n^nax = 6.12 n^at, M^nax = 2.263 Mq, R{M^^,) = 11.16 km] 
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c) NN force=SkI5 

[in npe matter in /3-equilibrium: UBic^ = 1) = 6.377 ng^t, n^lA Mq) = 2.10 Ugat, 
R{IA Mq)=13.87 km, n^^ax = 6.08 n,at, M^^, = 2.273 Mq, R{M^a^) = 11.36 km] 
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d) NN force=SV 

[in npe matter in /9-equilibrium: nsic^ = 1) = 4.983 ngat, ndlA Mq) = 2.10 ngat, 
R{IA Mq)=13A6 km, n^^ax = 5.44 n^at, M^nax = 2.426 Mq, R{M^a,) = 11.54 km] 
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Table V: Conditions for causality and neutron star properties 

In Fig. 6 we compare the mass-radius relation obained in several non relativistic models of 

matter in /3-cquilibrium with hypcrons. As discussed above, among our Skyrmc parametrizations 
the best arc the SkI3-hYBZ6-FSLL2 and the SLylO+LYI+SLL2 sets. The SkI5 or SV forces again 
give results very similar to those of the SkI3 force. The maximum mass reached with the stiffer 
SkI3-|-YBZ6-|-SLL2 set is 1.64 Mq comfortably above the Hulse-Taylor value 1.44 Mq while the 
maximum mass permitted by the softer SLylO-|-LYI-|-SLL2 set, 1.42 Mq, falls a bit short. This is 
however not uncommon a feature in non relativistic models (sec Brucckncr-Hartrcc-Fock calcula- 
tions [27, 28]). Among the parametrizations given by Balberg and Gal [2], that with intermediate 
compressibility corresponding to 7 = 5/3 was chosen. Figure 6 also displays the results obtained 
with a density-dependent Seyler-Blanchard potential as parametrized by Banik and Bandhyopadyay 
[24]. In contrast to our Skyrme forces the models of Balberg and Gal and of Banik and Bandhy- 
opadhyay both take into account S hyperons. The threshold for hyperon formation is represented 
on the figure with circles for the A and crosses for the S. 
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Fig. 6 - Neutron star mass-radius rela- 
tion for several non relativistic models of 
baryonic matter. 
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The parametrization of Balberg and Gal is in fact unstable slightly above the A hyperon thresh- 
old. As a consequence, the M vs. R relation displays an extended plateau, then recovers, crosses 
the line M = 1.4 Mq and finally reaches a maximum at M = 1.53 Mq. 

Also represented on the figure are lines of constant gravitational redshift 

^ 2GM\-^/2 ^ 



i?c2 ; 

To date two determinations exist from the observation of spectral lines in isolated neutron stars 
[29]. Sanwal et al. obtained a result with a large error bar z = 0.12 — 0.23 while Cottam et 
al. could extract a precise determination z = 0.35. All the neutron star models displayed would 
be compatible with the determination of Sanwal et al.. Almost all models are also in agreement 
with the value z = 0.35 of Cottam et al.. On the other hand, the mass-radius relation from the 
SkI3-|-YBZ6-|-SLL2 set is only marginally compatible with this value. 

Figure 7 displays the density profile of a 1.4 Af© neutron star and its hyperonic content. The pro- 
files corresponding to two of the Skyrmc parameters sets studied in this work, SkI3-|-YBZ6-|-SLL2 
and SLylO-|-LYI-|-SLL2 (top panels) are compared to a version of the models of Balberg and Gal 
and of Banik and Bandhyopadhyay (medium and bottom left panels) where all hyperons save the 
As arc artificially switched off. The profiles corresponding to the latter two models when also E~ 
hyperons are taken into account is displayed on the medium and bottom right panels. We can see 
that additional hyperons make the equation of state softer and the neutron star more compact. 
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Density profile of a 1.4 Mq neutron star for several non relativistic models of baryonic 



6 Finite temperature effects, neutrino trapping 

According to proto- neutron star formation and cooling calculations (see e.g. [30]), temperatures of 
the order of 30 to 50 MeV can be reached in the late phases of the supernova collapse. We therefore 
investigate here the range T G [0 — 50] MeV. 

The thermodynamical quantities are now to be written in terms of Fermi integrals. The expres- 
sions which are given in the Appendix at T = are easily generalized: The densities pi and the 
quantities Tj related to the kinetic energy should be replaced in the expressions for the baryonic 
contributions to the energy density and the effective masses by (in units h = c = ks = 1)' 



Pi = ^(2m*T)^/^/i/2(7?,) 

with In= 

Jo 1 + e«-^i 
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while the chemical potentials is related to the rji by fXi = rjiT + Ui{pi,Ti). The pressure is obtained 
from the derivative of the free energy 

T = £- — TS , ^ — ^2 + ^leptons , P = 1 2^"^ 

A,B=N,A P 

with the entropy 

S = Si + -Sieptons ) Si = Q^*rp ~ Pi Vi 

i=n,p,/V ' 

We used the GFD_D3 code published by Gong et al. [31] to calculate the Fermi integrals and 
keep as before the leptons fully relativistic whereas the baryons are treated nonrelativistically in 
consistency with the use of the Skyrme interaction. 



Skl3+YBZ6+SLL2 Skl3+YBZ6+SLL2 




Fig. 8 - Composition of neutron star matter in /3 equilibrium at finite temperature, calculated 
with the parametrization SkI3+YBZ6+SLL2 - (a) without neutrino trapping - (b) with neutrino 
trapping, Yl = 0.4 



The effect of temperature on the effective masses begins to be significant at T > 30 MeV; it 
tends to increase the neutron mass and decrease the proton effective mass so that their difference is 
reduced. The presence of trapped neutrinos, which tends to render the matter less asymmetric, also 
reduce the proton-neutron mass difference. In a warm protoneutron star, both effects are present 
and cumulate, the main contribution coming from Yi, ^ 0. The hyperon mass on the other hand is 
not significantly modified by temperature or neutrino trapping. 

The chemical composition is plotted on Fig. 8, without neutrino trapping on the left panel and 
with neutrino trapping on the right panel, for four values of the temperature, T=0, 10, 20 and 50 
MeV. We used in this section the parametrization SkI3 + YBZ6 + SLL2. 

Let us first discuss the neutrino-free case. It can be seen on Fig. 8-a that finite temperature 
effects are more important at moderate densities ub < 3nsat- Until T = 1 MeV the results are 
undistinguishable from the T=0 case. When T is increased from 1 to 50 MeV, the matter is more 
symmetric {i.e the proton fraction increases). There is stricto sensu no threshold for A hyperon 
production anymore, rather there always exist a vanishingly small number of As for any value of 
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the baryonic density below the T = threshold. Nevertheless, at the temperatures considered here, 
it is still possible to define a threshold density for practical purposes, which then moves to lower 
values as the temperature increases. 

The equation of state is somewhat stifFer at finite T, but this hardly affects the structure of the 
star until T=20 MeV. This result concern the case when no neutrinos are present in the matter. It 
is known however that high temperature shorten drastically the mean free path of the neutrinos so 
that a large amount of v are trapped. A typical value resulting from supernova collapse calculations 
is that the lepton fraction is of the order of Yl ~ 0.4. Fig. 8-b was drawn assuming this value for 
Yl. 

As expected the matter is more symmetric when neutrinos are trapped and the threshold 
for hyperon production is shifted to higher density. For example with the SkI3+YBZ6+SLL2 
parametrization and at T = we have (nthr = 2.08, Y*^'^ = 0.14) for = and (rithr = 2.84, 
Yp^^ = 0.36) for Yl = 0.4. For a given lepton fraction Yl, varying the temperature has even 
less impact on the composition of neutrino-trapped matter than in neutrino-free matter. Fig. 9-a 
shows that the neutrino trapped matter with hyperons is stiffer than its neutrino-free counterpart, 
a known result which leads to interesting consequences regarding a possible category of metastable 
stars which would collapse to a black hole as they cool when the deleptonization phase is completed 
(see e.g. [32]). Fig. 9-b illustrates this feature for the parameter set SkI3+YBZ6-|-SLL2: A newly 
formed protoneutron star with T = 30 MeV and Yl = 0.4 is metastable if its mass lies in the range 
M e [1.64 - 1.88] Mq (and radius in the range R G [12.08 - 14.07] km). For a lower mass, the star 
contracts as it looses its neutrinos and cools and its hyperonic content increases. 
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Fig. 9 - Effect of neutrino trapping and finite temperature on (a) the equation of state (b) the 
structure of the neutron star 
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As a conclusion to this subsection, our results concerning the effects of temperature and neutrino 
trapping are in full agreement with those obtained with other models of the baryonic interaction. 
The influence of the temperature on the equation of state and chemical equilibrium comes mostly 
indirectly through the buildup of an important fraction of trapped neutrinos. 
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7 Conclusion 



The aim of this work was threefold: 

- Use the most recent Skyrmc paramctrizations with hyperons existing on the market which are 
adjusted to reproduce the data on nuclei and hypernuclei and test them on neutron stars, 

- Study the influence of the hyperons on the ferromagnetic transition, 

- Inquire whether the hyperons are likely to affect the tail of the neutrino burst in super- 
nova explosions and prepare the background for the calculation of the neutrino mean free path in 
protoneutron stars. 

While the model considered in this work is still very schematic, it has led to several interesting 
results. 

The first result is that the presence of hyperons generally delay the onset of the ferromagnetic 
instability, and in many cases they are even able to remove it completely. 

An other advantage is that, by softening the equation of state, the hyperons remove the causality 
flaw and keep < 1 even up to very large densities. 

It is rather encouraging that NN and NA interactions from different authors and apparently 
very dissimilar parameter sets (compare e.g. SV to SkI3 or SKSHl to LY-I in Tables I, II !), once 
a series of reasonable requirements are fulfilled, yield very similar results qualitatively and even 
quantitatively. The major incognita is the A-A interaction which was very poorly known at the 
time the interactions used in this work were designed. 

After studying 43 NN forces in combination with 13 paramctrizations of the NA force and 4 
options for the AA forces, all taken from the litterature and known to reproduce correctly the prop- 
erties of nuclei and hypernuclei, four combinations of NN+NA+AA parameter sets were selected: 
(SkI3+YBZ6+no AA), (SkI3+YBZ6+SLL2), (SV+YBZ6+SLL2) and (SLylO+LYI+SLL2). Re- 
placing SkI3 by SkI5 or LYI by LYIV in the above sets would only give slightly different results. 
These sets fulfill all of the following conditions: 

(i) The NN force belongs to the subset selected by Rikovska Stone et al. The effective masses 
behave smoothly in all the relevant density and temperature range, pure neutron matter is always 
stable, the asymmetry energy does not decrease so much with density that protons would disappear 
from the system, neutron stars formed from npe matter in P equilibrium can reach a mass at least 
equal to 1.4 Mq 

(a) The NA force belongs to the set preferred by Lanskoy et al. as best reproducing the 
properties of hypernuclei. The threshold for hyperon formation should lay between 1.7 and 4 times 
saturation 

(in) The neutron star formed of softer npAe matter in (3 equilibrium should still reach a mass 
at least equal to 1.4 Mq. 

(iv)'No ferromagnetic transition should be present As explained in the main text, this is, for a 
Skyrme parametrization, a strong requirement. 

The second result of this work is that the selected sets are also found to reproduce all features 
observed in other models such as non-relativistic Brueckner-Hartree-Fock or relativistic mean field 
calculations, not only qualitatively but also quantitatively. The threshold for hyperon formation is 
in fact restricted to the narrower range [2-2.5] ngat- The softening of the equation of state brings 
the maximum mass of the star from 2 - 2.5 Mq for a npe star down to 1.4 - 1.6 AIq for a npYe 
star. This has even lead to some speculation (see e.g. [27]) whether the clustering of known pulsar 
masses around the value 1.4 Mq would be due to this feature of the equation of state rather than 
from the circumstances of their formation in a supernova. The metastability of hot stars with 
trapped neutrinos when hyperons are present, as discussed e.g. in [32], is also recovered. 

The selected sets are also applied in a companion paper [16] to the calculation of the scattering 
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rate of neutrinos in baryonic matter, for a finite temperature and non vanishing number of trapped 
neutrinos. The possibihty to access to the properties of polarized matter was a further motivation 
to this work: It was necessary to be able to calculate the Landau parameters in the spin S = 1 
channel, since the axial channel is dominant for neutrino scattering. 

As a protoneutron star is formed in a supernova explosion, the high amount of neutrinos trapped 
in its interior by high temperature and density gradually diffuse out in the first 50 seconds as the 
star cools. Hyperons begin to appear in the course of this deleptonization process and could affect 
the tail of the neutrino signal. Modern detectors would now be able to detect this effect. The idea 
was explored by Reddy et al [33] in the mean field approximation. The calculation is performed in 
the random phase approximation in [16]. The parameter sets selected in the present work predict 
a non negligible hypcron content for t > 20 s after the collapse (sec Figs 2 and 3 of [16]). 

Let us examine again the main shortcomings of the present model and the way it could be 
improved in future work. 

(i) The Skyrme interaction is long known for displaying a series of problems at high density: 
behavior of asymmetry energy, onset of fcrromagnetism, causality ... While we were able to avoid 
these problems by carefully passing the existing parametrizations through the crib of these various 
constraints, we cannot avoid the feeling that we are pushing the Skyrme model much beyond 
its capacity. More refined models do not have these problems; for example the non relativistic 
Brueckner Hartree Fock calculations do not show a ferromagnetic transition. A possible answer 
to this point would be to develop a parametrization of Brueckner-Hartree-Fock calculations in 
terms of an energy density functional for polarized matter with a non vanishing hyperonic content. 
Some steps have already been performed in this direction by Vidaha et al. These authors studied 
the polarized neutron-proton matter system [12] and concluded to the absence of a ferromagnetic 
transition. They also obtained the equation of state of unpolarized npY matter in (3 equilibrium and 
applied it to the calculation of neutron star mass radius relation [27] and parametrized the neutron- 
proton-Lambda system [34] for applications to hypernuclei. The parametrization of Brueckner- 
Hartree-Fock results with modern Nijmegen potentials for the full polarized npAS are currently 
under way [35]. 

(ii) We have seen that the effect of the A- A force at high density is very important, whereas the 
description of the hyperon-hyperon interaction in phenomenological formalisms is still very poor. 
Brueckner-Hartree-Fock calculations offer a coherent framework to calculate the hyperon-hyperon 
interaction in medium starting from bare potentials such as the Nijmegen one known from scattering 
data of free particles, and then cross checking the results with the data on double hypernuclei. This 
has actually been performed with the Nijmegen NSC97e model e.g. in [36]. 

(in) Relativistic effects can be expected to play an important role at high density. Many of the 
problems encountered in the Skyrme parametrization mentioned above in this paragraph (behavior 
of asymmetry energy, onset of ferromagnetism and causality) are not present in the relativistic 
formulation. Relativistic mean field models with hyperons (see e.g. [37]) are so well under control 
that they have made their way into textbooks [38]. The parametrizations of the present work 
obviously do not pretend to compete with this line of work at the mean field level; rather they were 
developped in view of their application at the RPA level where they give rise to simpler Dyson 
equations than in the relativistic formulation. The relativistic extension at RPA level to hyperons 
would be completely straightforward but somewhat unwieldy. 

While it could be argued that the physical foundations of the Skyrme model are disputable 
when applied in the context of neutron stars, the parametrizations presented in this paper should 
rather be considered as reliable phenomenological models for run-of-the-mill calculations. Once the 
order of magnitude of physical effects are ascertained with this simple model, they can refined by 
applying more complicated Brueckner-Hartree-Fock or/and relativistic models. 
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8 Appendix 

Besides the the usual parameterization of the nucleon-nucleon interaction 

VNN{ri - ra) = to (1 + XQP^)5{ri - rs) + \h (1 + x^P^) \k'^ 5{n - r2) + 5{rx - ra) A; 



+ t2 (1 + X2Pa)k' 5{ri -r2)k + (1 + x^P^) p% 



(^^)5(ri-r2), (14) 



we use the following Lambda-nucleon and Lambda-Lambda potentials [4, 5] : 

VNA{rN - ^a) = uo{l + yoPa)6{rN - ta) + ^ui ^k'"^ 6{rN - va) + 6{rN - ta) fe^ 



+ U2 k' S{rN -rA)k + (1 + y3Pa)/4 (^^^^4^) ^^'''^ ~ ""^^ ^^^^ 



VaaIh - r2) 



1 



= Ao5(ri - r2) + -Ai l^/c 6{ri - r2) + S{ri - r2) k" 
+ X2 k' 5{ri -r2)k + A3 paPat 



(16) 



(We haved dropped in these expressions the spin-orbit terms which are not used in this paper) . We 
obtain the energy density in homogeneous matter in the usual way. In spin saturated npA matter, 
the functional reads: 



S 

£nn 



< iIj\H [iIj >= Snn + £na + ^AA 
2mi\f 2 



xo 



1 



1 + ^ - o + (Pn + p: 



l + ^)p'N-{^+Xs\{pi+p^) 



, *3 a 

+ 1 [ti(2 + Xi) + t2(2 + X2)] PNTN - I [tl{2xi + 1) - t2(2X2 + 1)] (pn^n + PpTp) 
o o 



(17) 



^ATA 



^AA = 



yo 



^^0 (^1 + yj PivPA + g'ua Pat^Va + ^ 
+^ [ui{2 + yi) + ti2(2 + y2)] (patta + Patat) 



2m\ 



+ yPA + ^(^1 + 3A2) PATA + ^ PAPAf 



A3 
4 



(18) 
(19) 



with pn = Pn + Ppi tn = Tn + Tp. At T = we have ti = (3/5) Pikp^^ and pi = kp^/ (Stt^). At T 7^ 
these expressions should be replaced by Fermi integrals, see §6. At T = we obtain the chemical 
potentials 



Pn 

Un{pi,Ti) 



ds _ 

dpn 2m. 

to 



+ 24PN 



[p„(l - Xo) + Pp(2 + Xo)] 
pI {(2 + a)(l - X3)} + pI {2 (2 + X3) + a{l - X3)} + 2pnPp {3 + a(2 + X3)} 



+ ^ [ti(l - xi) + 3t2(l + X2)] r„ + ^ [h{2 + xi) + t2(2 + X2)] Tp + ^-/pIpI ^ 
+ «o (1 + y) PA + lu3{(5 + 1) (1 + ^) p^PA + I [ui{2 + 2/1) + «2(2 + ^2)] ta 



(20) 
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(The chemical potential for the proton /Xp can be obtained from fj,^ by interchanging the indices n 
and p) and 

dS ii , , / N 

UtXpi.Ti) = yPa + Yp]irPA + ^(Ai + 3A2)rA 

+ uo(l + ^) PN + Im (l + ^) p^' + \[ur{2 + yi) + U2{2 + y2)]TN (21) 



2 ; " 8 V 2 
The effective masses are given by 



+ - [ii(2 + Xi) + i2(2 + X2)\ PN-^ [il(2xi + 1) - t2{2X2 + 1)] p„ 



2m* 2mAr 8 

+\[ui{2 + yi) + U2{2 + y2)]pK (22) 

(The effective mass of the proton follows from replacing p„ by pp in this expression) and 

1 1 
+ - [Ai + 3A2] PA + - [ui{2 + yi) + U2(2 + 2/2)] PN (23) 



2m^ 2mA 8 

For deriving the Landau parameters we will also need the polarized energy functional £{pni,Pni, 
Pph Ppi^ P^h P^i)- reproduce this somewhat lengthy expression here, but note that it 

coincides for np matter with the functional given by Bender et al. [39]. 

The Landau parameters in the monopolar approximation ^ = can be calculated in the following 
way. We first taJie the derivative of the single particle energies (or equivalently the second functional 
derivative of the energy density with respect to occupation numbers piik)): 

(24) 

In the Skyrme model the single particle energies U-j-aik) are quadratic in the momentum k 

Ura{k) := Urc + ^k^ (25) 
and related to the chemical potentials in polarized matter 

l^'^'^ •= = + —ir^FTa = Uraikpra) (26) 

\(^PTa J \p{T'<T'i^Ta)=CSt '^TC 

with pra = kl^^/iGTT^) at T = 0. 

The Landau parameters in the spin S = channel are obtained from 

/rir2 ■= ^ [/nTT-2T + /nj.r2i + /TiTT2i + /TiiT2T]|unpolarized (^'^) 

and in the spin S = 1 channel 

9tiT2 T [/TitT2T "I" /tiJ.i"2J. ~ /TitT2j. ~ iT-2 tI |unpolaxized 

(28) 
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The fj-^j-^ are related to usual thermodynamical quantities, for example the compressibility and the 
asymmetry energy 



K 
with 



m 



*asym 



N 



6m 



N 



kp — — kpp , fnn — fpp j /o 



fpp + /) 



np 



) /o 



fpp fi 



np 



Fo = No /o, F^ = iVo /^,A^o 
in symmetric nuclear matter. 



^2 



(29) 



The gTiT2 are related to the magnetic susceptibilities and are used to form the ferromagnetic criterion 
(Eqs. 11 - 13) 

Their explicit expressions are 
In the spin S = channel 

1 



fn 



to{l - xq) + -^tspffil - X3) + -ahp% ^ 



12 



(1 + ^)PN 



(2 + ^3)Pr, 



1 



2 



(1 + vK-(o+^3)(Pn + Pp) 



Inp 



+^n3(l + |)/3(^ + 1)p^-Va + ^A37(7 - l)prVi 
io(l + y ) + ^t3P^(l + y ) + ^atsp^ 



+ -[ti(l-a;i) + 3*2(1 + 0:2)] 4n 

(30) 



+-a(a - l)t3Pr' 



{kin + kip) 



1 

+ 4 



ii(l + y)+^2(l + f 



+^n3(l + |)/3(^ + 1)p^-Va + ^A37(7 - l)pli'pl 
^to{l - xo) + ^tspffil - X3) + ^at3p% 



(1 + y )Piv -{\ + x?,)pp 



(31) 



(l + yK-(^+:^^3)(p^ + /^^) 



1 



-[ti(l-xi) + 3*2(1 + 0:2)] kip 



+^n3(l + |)/3(/3 + 1)p^-Va + ^A37(7 - l)p]rVi 



/nA — 



13 1 

-uo(2 + yo) + Y^^3(2 + y3)(l + /3)P^ + g [^^i(2 + yi) + U2{2 + 2/2)] (4a + 4„) 



1 



+^A37p]^ Va 



fpK = -UQ{2 + yQ) + —U3{2 + y3){l + P)p% + -[ui{2 + y^)+U2{2 + y2)]{kl^ + klp) 



/aa = 



2 16 

+2^3^p/rVA 

^Ao + \\3pl + J [Ai + 3A2] 4a 



(32) 

(33) 

(34) 
(35) 



In the spin S = 1 channel 



^to(aJo - 1) + ^hp%ix3 - 1) + J [hixi - 1) + i2(l + X2)] 4„ 



(36) 
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gnp = ^toXo + —t3X3p% + - [tiXi + t2X2] {kp^ + kpp) (37) 

9pp = \to{xo-l) + ^t3P%{x3-l) + ^[h{xi-l) + t2{l + X2)]kj^p (38) 
13 1 

gnA = ^uoyo + Yg^sysPAT + g [""lyi + ^^2^2] {kpj^ + A;p„) (39) 
1 3 1 

9pA = ^uoVo + Y^uayspfj + - [myi + ^2^2] (A^fa + ^Fp) (40) 

9AA = ~>^o-^X3pJf + \[-^i + ^2]klj^ (41) 

In the limit where hyperons are absent these expressions coincide with those of Hernandez et 
al. [40] 
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